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Subject of this letter is the dynamics of a chain obtained performing the continuous 
hmit of a system of hnks and beads. In particular, the probability distribution of 
the relative position between two points of the chain averaged over a given interval 
of time is computed. The physical meaning of the obtained result is investigated in 
the limiting case of a stiff chain. 



INTRODUCTION 



We study in this letter the dynamics of a chain obtained performing the continuous limit 
of a system of links and beads. This kind of problems has been addressed in the seminal paper 
of Edwards and Goodyear l| using an approach based on the Langevin equation. Related 
papers dedicated to the statistical mechanics of a freely jointed chain in the continuous limit 



are for example 



mm. 



In Ref. 



6j it has been shown that it is possible to investigate 
the dynamics of such a chain in a path integral framework. The model which describes 
the fluctuations of the chain is a generalization of the nonlinear sigma model [3] called the 
geuera.i.ed nouMe.. sigma mode, or simply GNLSM. I„ H it has also been diseussed the 
relation between the GNLSM and the Rouse model [8|, |9|. Applications of the GNLSM have 
been the subject of Ref. [l^, in which the dynamic structure factor of the chain has been 
computed in the semiclassical approximation. 

The goal of the present work is to derive the probability distribution Z{ti2) which mea- 
sures the probability that in a given interval of time At the average distance between two 



2 



points of the chain Pi and P2 is ri2. With respect to Ref. [10(], we use to perform the calcula- 
tion a totally different approximation, which linearizes the GNLSM. Moreover, a background 
field method is adopted, in which the effects of the thermal fluctuations are considered in 
the background of a fixed chain configuration chosen among the classical solutions of the 
linearized equations of motion. We find out that classical configurations are particularly 
stable against the changes due to fluctuations. The latter become relevant only if they act 
on the chain for a significant amount of time. 

The material presented in this letter is divided as follows. The GNLSM is linearized 
exploiting a gaussian approximation of the functional Dirac delta function. Next, the prob- 
ability distribution Z{ri2) is computed within this approximation. The asymptotic form of 
Z{ri2), which is valid in the case in which the chain is stiff, is derived. Finally, a discussion 
of the obtained result is presented. 

THE GENERALIZED NONLINEAR SIGMA MODEL IN THE GAUSSIAN 

APPROXIMATION 

Let us consider the partition function of the GNLSM: 

Z = J PR(t,(T)e-'^/o'r'='-^'(*'-)5(|R'| -£) (1) 

with R = dR/dt and R' = dH/da. The boundary conditions at t = and t = tf oi the field 
R are respectively given by R(0, a) = Ro(cr) and R(t/, cr) = R/(cr), where Ro(o") and R/(cr) 
represent given static conformations of the chain. The boundary conditions with respect to 



a are periodic: R(t, cr) = R(t, a + A^). It was shown in Refs. p] and lOj] that the above 
partition function describes the dynamics of a closed chain that is the continuous version of 
a freely jointed chain consisting of links and beads. The constant c appearing in Eq. ([1]) is 
given by: 

c = ci with c = — (2) 

Here denotes the Boltzmann constant, T is the temperature and r is the relaxation 
time which characterizes the ratio of the decay of the drift velocity of the beads. M and L 
represent the total mass and the total length of the chain respectively. 



Let us note that in Eq. ([T]) the trajectory of the chain has been parametrized with t 
help of the dimensionless parameter a, which is related to the arc-length used in Refs. 



le 



and 1^ by the relation s = ia. The scale of length i introduced with this parametrization 
is connected to the quantities L and by the identity N = j. We remark also that, with 
respect to Refs. jd] and 10 1, the constraint R'^ = £^ has been replaced by the constraint 
|R'| = i exploiting the property of the functional Dirac delta function 5(R'^ — = 5(|R'| — 
£). A proof of this equation can be found in Ref. Q]. 

In order to deal with the delta function in Eq. ([1]) we use the following gaussian approx- 
imation [gI: 

5i\R'\ -i)r^ exp f /*' dt r rfa^R'2 ) (3) 



which is valid when the parameter v is large, while £ is small. As a consequence, the partition 
function Z becomes: 



PR exp 



N 



dt I da (gR2 + -R'2 
JO V 2 



(4) 



After the approximation ([3]) the chain may be regarded as a gaussian chain consisting in a 
set of A^ segments of average length i. 



THE PROBABILITY DISTRIBUTION OF THE AVERAGE POSITION 
BETWEEN TWO POINTS OF THE CHAIN 



At this point we pick up two points Pi and P2 of the chain, for instance R(t, ai) and 
R(t, 0-2). We wish to study how the relative position between these two points changes due 
to thermal fluctuations. To this purpose, we introduce the following probability distribution: 



Z(r 



12 



I^Rexp 



dt r da f5R2 + ^R'2 
JO V 2 



X 



'5(ri2-££(R(t,a2)-R(t,ai))) 



(5) 



We suppose that < ti < ^2 < To evaluate the path integral in Eq. ([5]) it is convenient 
to perform the splitting: R(t, a) = Rci(t, cr) + rit.a). Here Rc/(t,cr) is a solution of the 



classical equations of motion + |^ Rd = 0. The fluctuations around the classical 
background Rd(t, a) are taken into account by r(t, a). The boundary conditions for Rc/(t, a) 
at the initial and final instants are: Rcz(t/, a) = R/(o") and Rd(0, a) = Ro(cr). With respect 
to a periodic boundary conditions are assumed: Hciitj, a) = Hd(t, a + N). The solution of 
the classical equations of motion is complicated due the presence of the non-trivial boundary 
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conditions and will not be reported here. It can be found in standard books of mathematical 

n 

methods in physics, like for instance [11]. The only important information related to the 
background fields that is needed in the probability distribution Z{yci) is the average relative 
position Yd of the points Pi and P2 for a given classical conformation Rd(t, cr): 

r*2 dt 



rt2 

rci = ^ (R-ci{t, 0-2) - Rci{t, 0-1)) 



(6) 



For the fluctuations r(t,cr) we choose Dirichlet boundary conditions in time r{tf,0) = 0, 
r(0, cr) = 0. The boundary conditions with respect to a are instead periodic: r(t, a) = 
r(t,a + N). After a few calculations, it is possible to rewrite Z{vi2) in the form: 



/oo 
-co 



(7) 



where r^i and = Jq^ dt da (cTlli + |Rcf ) take into account the contribution of the 
classical background, while the fluctuations appear in Z{k): 



Z(k) = Jvvexp - J^^ dt da (^cr"^ + '^r''^ + iJ ■ r 



(8) 



In Eq. dH]) the external current J is given by: 



Jit, ^) = ^ (^(^ - ^2) - S{a - a^)) e{t2 - t)e(t - ti) 



(9) 



This current has been introduced to provide a convenient way of rewriting the integral 
/t^f(r(t,^2)-r(t,ai)). 

Let's now concentrate on the computation of the partition function of the fluctuations 
Z(k). The gaussian integration over the flelds r(t, a) may be easily carried out and gives as 
a result: 

Z(k) = Ce^C') (10) 



where C is a constant and 



VF(k) = exp 



- / dtdt' / dada'G{t, a; t', a')3{t, a) ■ J(t', a') 
2 JO Jo 

In the above equation G{t, a; t', a') denotes the Green function satisfying the equation: 

G{t,a;t',a') = 6{t-t')6{a-a') (12) 



(11) 
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According to our settings, we choose for G{t, a; t', a') Dirichlet boundary conditions at the 
instants t = tj and t = and periodic boundary conditions in a and a'. To solve Eq. ([T2l) . 
we decompose G{t, a; t', a') in Fourier series as follows: 

+00 

Git,a;t',a')= 9nMe-^^^i^^~^'> (13) 

n=— 00 

Substituting also the Fourier expansion of the periodic Dirac delta function 5{a — a') = 
En=-oo ^e-'^'^'T in Eq. ^ and solving for gn{t,t'), we obtain for n = 0: 

9o{t, t') = ^e{t' - t)t^ + ^e{t - t')t'^—^ (14) 

LiC if LC If 



and for n 7^ 0: 



+ A Jit' - t) sinh I3nt sinh /?„(t/ - f) (15) 



with 



= (16) 
= 7^ , , .\ , (17) 

V2z/^c |n|7rsinh/^„t/ 

In Eqs. (JUD and <^ 0{t) denotes the Heaviside 6'-function 0{t) = 1 for t > and 9{t) = 
for t < 0. To complete the calculation of Z(k) in Eq. (fTOj) we rewrite W^k) as follows: 

1 +00 „f „jy 

W{k) = - y / dtdt' d(jda'gn{t,t')e-^^'Li---')^:i{t,a)-:i{t\a') (18) 

It is easy to show that the only non-zero contributions to the above integral are those for 
which n 7^ 0. After integrating over the variables cr, a' and over one of the time variables, 
we obtain: 



wik) = y — 



1 — COS 27m 



'2 



L 



X 



ri2 

/ (it sinh /3„(i(:/ — t) (cosh /9„t — cosh /?„ti) (19) 

J ti 



where we have put £2 = ^cr2 and ii = iai. Remembering that the coefficients An defined 
in Eq. ( JT71) are all strictly negative, it is easy to realize that iy(k) is negative too. This 
property of W^(k) will be necessary in order to perform the remaining integration over k in 
the probability distribution Z{ri2) of Eq. ([7j). 
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A last integration over dt in Eq. (fT9|) delivers the following result: 

'1 



iy(k) 



1 — cos 27m 



/ X . , ^ cosh. But f 



4/?n 



(cosli/3„(t/ - 2t2) + cosh/5„(t/ - 2ti)) 



(20) 



+ — (cOsll/?„(tj+ti-t2)+COsh/?„(tj- (ti+t2))) 



Substituting Eq. fl20l) back in Eq. ffTOj) we get an exact expression for 2'(k). To obtain the 
probability distribution Z^r^) of Eq. (JTj) in closed form a gaussian integration over k is 
sufficient. 



CALCULATION OF THE PROBABILITY DISTRIBUTION ^(ria) IN THE LIMIT 

OF A STIFF CHAIN 

To simplify our calculations of the previous Section, we compute the infinite sum in 
Eq. (!20|) using an approximation. First of all, we give a physical meaning to the parameter 
u by putting: 

After performing the above substitution in the partition function of Eq. (jlj), the term pro- 
portional to u appearing in the action of the fields R becomes exactly the term introduced 
in the GNLSM in Ref. in order to take into account the bending energy of the chain. 
This connection with the bending energy is confirmed by the fact that a has dimensions of 
an energy per unit of length, i. e. [a] = [energy] ■ [length]"^. 



As it was shown in Ref. [10(], large values of a correspond to a stiff chain. Indeed, it is 
possible to check from Eq. ([8]) that in this case the corrections to the tangent vectors R'^; 
coming from the fluctuations r' are strongly suppressed in 2'(k). In the following, we will 
work in the limit of a stiff chain, i. e.: 



a » 1 (22) 

implies that also From Eqs. (!2T!) and (!22l) it turns out that the quantity ui = appearing 
in the coefficients /5„ and An (see Eqs. f|T6l) and f|T7|) ) is very large. As a consequence, 
within the approximation (122!) the coefficients are large, while the coefficients An are 
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exponentially small. Taking into account these facts, we may write the following asymptotic 
expression of W^(k): 

2k2 L 



(At)2 27r2i 



c 

TT V 2^ 



(23) 



In Eq. 



we have introduced the convenient notation: 



B 



1 



E 



1 — cos 27m 



1 — cos27rn 



L 



(24) 
(25) 



It is possible to check that the other contributions to Vr(k), expressed in Eq. (l23l) with the 
symbol O (^~~^f^^*'^^^^^^^ decay at least as fast as e~~v^~^''*^~*^\ These terms become 
negligibly small when vl is large, provided of course that: 



IL re 



t2-tl » -- 

/ TT V lui 

Substituting Eq. (I23l) in Eq. ffTOj) . we obtain for 2'(k) the approximate expression: 



(26) 



Z{k) = Ce 



2 



where 



K 



(At)2 7r22z/£ 



B 



L 



(27) 



(28) 



2 7r V 2ui\ 

The condition (l26l) guarantees that k, > 0. This is what is needed to compute the integral 
over k in Eq. ([7j). After a few calculations we arrive at the final result: 

1 (ri2 - r^if' 



Z{ru) = Ce-^- ( ^) ' exp 



K 



(29) 



Eq. fl29|) has a straightforward interpretation. The quantity k is inversely proportional to the 
product ui, which is supposed to be large. Thus k is very small. As a consequence, Eq. (!29|) 
implies that the relative positions of the two points Pi and P2 exhibits a sharp peak around 
the classical value, i. e. when ri2 = r^/. From Eq. fl29p it turns also out that the changes of 
the classical background configuration Rc/(t,cT) due to the fluctuations r(t, a) are relatively 
small. This is due to the fact that there is no contribution to the parameter k, which is 
of zeroth order in ui. Potentially, a zeroth order contribution could come from the term 
proportional to go{t, t') in Eq. (fT8|) . but we have seen that this term vanishes identically. This 
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suggests that the effects of thermal fluctuations are weak at short time-scales. Of course, 
with the increasing of the time interval At = t2 — ti over which the average of the relative 
position of Pi and P2 is measured, the fluctuations have more and more time to act on the 
chain and their influence becomes more and more important. This intuitive prediction is 
confirmed by Eq. ( |29l) . Indeed, the coefficient k, grows linearly with At and, for larger values 
of K, the peak around the classical point ri2 = r^i in the probability distribution Z{yi2) 
becomes less sharp as expected. 



CONCLUSIONS 



In this work the GNLSM of Refs. M and 



.a 



lOj has been applied to compute the distribution 



function Z{yi2) of the relative position between two points of a chain with rigid constraints. 
The calculation has been performed using the approximation ([3]) of the functional Dirac 
delta function which is needed to impose the constraints. After this approximation, the 
finest details of the chain are lost, a fact that has already been noted in the case of a static 
chain 

A closed form of the probability distribution fl29l) may be obtained starting from the 
exact calculation of the contribution of the fluctuations to Z{yi2) given in Eq. (!20l) . The 
final formula for Z{ri2) computed in this way is however complicated. For this reason the 
physical meaning of Z{ri2) has been investigated at the end of the previous Section in the 
case in which the energy needed for bending the chain is large. In this approximation 
the contributions of the fluctuations which decay exponentially when a becomes large are 
neglected. 

Concluding, it would be interesting to explore the connections between the GNLSM and 
other models of the dynamics of a chain, in which instead of rigid constraints the addition 

n 

of potentials is used in order to prevent the breaking of the chain [13l |. 
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